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Abstract. We produce lensing potential and deflect ion- angle maps in order to simulate the weak 
gravitational lensing of the Cosmic Microwave Background (CMB) via ray-tracing through the 
COupled Dark Energy Cosmological Simulations (CoDECS), the largest suite of N-body simulations 
to date for interacting Dark Energy cosmologies. The constructed maps faithfully reflect the N- 
body cosmic structures on a range of scales going from the arcminute to the degree scale, limited 
only by the resolution and extension of the simulations. We investigate the variation of the lensing 
pattern due to the underlying Dark Energy (DE) dynamics, characterised by different background 
and perturbation behaviours as a consequence of the interaction between the DE field and Cold 
Dark Matter (CDM). In particular, we study in detail the results from three cosmological models 
differing in the background and perturbations evolution at the epoch in which the lensing cross 
section is most effective, corresponding to a redshift of ^ 1, with the purpose to isolate their 
imprints in the lensing observables, regardless of the compatibility of these models with present 
constraints. The scenarios investigated here include a reference ACDM cosmology, a standard 
coupled DE (cDE) scenario, and a "bouncing" cDE scenario. For the standard cDE scenario, we 
find that typical differences in the lensing potential result from two effects: the enhanced growth 
of linear CDM density fluctuations with respect to the ACDM case, and the modified nonlinear 
dynamics of collapsed structures induced by the DE-CDM interaction. As a consequence, CMB 
lensing highlights the DE impact in the cosmological expansion, even in the degenerate case where 
the amplitude of the linear matter density perturbations, parametrised through as, is the same in 
both the standard cDE and ACDM cosmologies. For the "bouncing" scenario, we find that the two 
opposite behaviours of the lens density contrast and of the matter abundance lead to a counter- 
intuitive effect, making the power of the lensing signal in this model lower by 10% than in the 
ACDM scenario. Moreover, we compare the behaviour of CDM and baryons in CoDECS separately, 
in order to isolate effects coming from the coupling with the DE component. We find that, in the 
bouncing scenario, baryons show an opposite trend with respect to CDM, due to the coupling of 
the latter with the DE component. These results confirm the relevance of CMB lensing as a probe 
for DE at the early stages of cosmic acceleration, and demonstrate the reliability of N-body based 
large scale CMB lensing simulations in the context of DE studies. 
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1 Introduction 

Anisotropics in the Cosmic Microwave Background (CMB) represent a pillar of our present under- 
standing of cosmology. They are sourced by all kinds of cosmological perturbations, namely scalar 
(such as density), vector (vorticity) and tensor (Gravitational Waves, GWs). At recombination, 
Thomson scattering stores imprints of such perturbations in the intensity and linear polarisation of 
photons and enables us to observe them as primary anisotropics in the CMB sky, to be distinguished 
from the secondary ones collected along the line of sight. The scattering produces total intensity 
(T) anisotropics and linear polarisation ones, described by the Stokes parameters Q and (7, see e.g. 
Hu et al. [1] for reviews. The latter are further combined into gradient modes (or electric-type, E) 
and curl ones (or magnetic- type, B). At the linear level, while the E modes are sourced by all kinds 
of perturbations, B modes are excited only by vorticity and GWs [2, 3]. 

The T power on large scales, greatly dominated by the metric effect known as Sachs- Wolfe (SW), 
occupies angular scales up to the angle subtended by the horizon at recombination, corresponding 
to a multipole / 2^ 200 in spherical harmonics; degree and sub-degree angular scales are domi- 
nated by acoustic oscillations from radiative pressure opposing gravity. Polarisation is caused by 
local anisotropy in T on the last scattering region, and has significant power on the degree and 
sub-degree scales, featuring acoustic oscillations, largely correlated with those in T. 

The Planck satellite [4] has observed CMB anisotropics for T over the whole sky down to a 
resolution of 5 arcminutes, with a sensitivity of a few /iK per resolution element, observing large 
scales and seven acoustic oscillations. Measurements of the TE correlation have also been reported 
[4] and the analysis of data is still ongoing, in particular the polarisation measurements will be 
made public in a forthcoming release. Previously, the Wilkinson Microwave Anisotropy Probe 
(WMAP) reached an angular resolution of about 20 arcminutes, also producing measurements of 
the TE correlation [see e.g. 5, and references therein]. Operating experiments from the ground and 
monitoring percent sky fractions have already reported the existence of acoustic oscillations for T 
down to the arcminute scale, namely the South Pole Telescope [SPT, see 6, 7] and the Atacama 
Cosmology Telescope [APT, see 8, 9]. Sub-orbital probes such as the E and B experiment^ [EBEX, 
see 10] and POLARBEAR ^ [11], Spider [12], SPTpol [13], ACTpol [14], the Large Scale Polarisation 
Experiment [LSPE, see 15], QUbic [16] are designed to probe CMB polarisation looking in particular 
for the B modes from primordial GWs on percent sky fractions with angular resolution similar to 
Planck; moreover, proposals for a future post-Planck CMB satellite, aiming at an all sky study of 
CMB polarisation, are in preparation: see e.g. the Cosmic ORigin Explorer^ (CORE), CMBpol^, 

^ groups . physics . umn . edu/ cosmology /ebex 
^bolo. berkeley.edu/polarbear 
^ www. core-mission.org 
^cmbpol. uchicago.edu 



the Primordial Inflation Explorer [17], and the Light (Lite) satellite for the studies of B-mode 
polarisation and Inflation from cosmic background Radiation Detection (LiteBIRD^). 

The data available at present are well represented by a six parameter cosmological model 
known as ACDM [see 18, and references therein]: the expansion proceeds with a rate of ^ 68 
km/s/Mpc, and is accelerating under the effect of a Dark Energy (DE) component constituting 
^ 68% of the entire cosmic energy density budget, the rest being described as Cold Dark Matter 
(CDM) responsible for the dark halos around galaxies and galaxy clusters (about 27%) as well 
as leptons and baryons (about 5%); the primordial density power spectrum departs from scale 
invariance, oc A:^^, where Ug ^ 0.96. No detection of cosmological GWs exists to date. 

Present and future CMB observations are constraining the cosmological dark components, 
through the interaction between the CMB and the forming cosmological Large Scale Structures 
(LSS) along the hne of sight. This latter aspect represents the context of the present work. 

Secondary anisotropics are caused by the interaction of CMB photons with LSS along their 
paths from the last scattering surface to the observer. CMB photons re-scatter onto electrons 
freed by cosmological reionisation at the beginning of structure formation. This boosts polarisation 
anisotropics on the angular scales subtended by the horizon at that epoch, corresponding to several 
degrees in the sky, as it has been detected by the WMAP satellite [see e.g. 5, and references 
therein]. Moreover, the Integrated SW (ISW) effect is caused by CMB photons crossing regions 
characterised by evolving gravitational potentials, and has been now observed by Planck [19]; this 
effect is active already for linearly evolving LSS while it is known as Rees-Sciama (RS) when the 
underlying structures are non-linear. The efficiency of the ISW reconstruction from data is boosted 
when CMB and LSS are Cross Correlated (XC); several authors have used techniques aiming at 
estimating such a correlation in order to extract cosmological constraints on DE as well as on the 
statistics of primordial perturbations from the ISW [20, 21, 22, 23]. 

The subject of this work concerns a third aspect of secondary anisotropics, namely the capa- 
bility of LSS to act as gravitational lenses for CMB photons. The combination of the statistical 
spatial distribution of power in the CMB anisotropics and LSS, associated with the geometry of the 
expansion, makes the CMB lensing process effective on sub-degree and arcminute angular scales, 
causing a smearing of acoustic peaks and the transfer of power in the damping tail of CMB primary 
anisotropics, as well as a transfer of the E mode power into the B modes, resulting in a characteristic 
peak in the latter at multipoles of about / ^ 1000 [24]. Planck has measured the CMB lensing with 
about 26cr confidence level, covering the whole spectrum of scales where it is effective [25], along 
its correlation with the structures observed as the Cosmic Infrared Background [CIB, see 26, and 
references therein]. Moreover, previously and independently, the ACT [9] and SPT [7] experiments 
detected the lensing signal in the damping tail of CMB total intensity anisotropics. 

The lensing cross section peaks half way between the observer and the source. If the source 
is effectively at infinity as in the case of CMB, the epoch in which the lensing distortion is most 
recorded into CMB anisotropics coincides with the onset of cosmic acceleration, at z 2^ 1 ± 0.5: 
therefore, since the lensing efficiency is related to the behaviour of the expansion and perturbations 
at the corresponding epoch, the CMB lensing signal carries important information about the onset 
of late-time cosmic acceleration, poorly constrained at present, promoting CMB alone to be a probe 
of the DE behaviour at redshift of about 1, regardless of how close it is to a Cosmological Constant 
at present [27]. On the contrary, primary CMB anisotropics suffer the projection degeneracy making 
them unable to distinguish a global curvature from a Cosmological Constant. Predictions say that 
CMB lensing should be able to measure the DE abundance at the onset of acceleration with ^ 10% 
precision [27, 28]. ACT and and SPT data have already been used [29, 30] in combination with 
other cosmological data sets, to break the projection degeneracy and give evidence of DE from the 
CMB alone. Recent theoretical works are investigating the capabilities of CMB lensing to measure 
not only the abundance of DE at the onset of acceleration, but its coupling with CDM and of early 
DE [31]. The CMB lensing as a DE probe will benefit from the XC with the observations of the 
actual lenses in LSS surveys. The XC with CMB lensing has been already detected on the basis 

^http:/ /cmbpol. kek.jp/litebird 
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of the available CMB and LSS data [32, 33]. The latter will culminate in about a decade with the 
observations of the Euclid satellite^, which will perform arcsecond imaging of billions of galaxies 
over half of the sky between redshifts and 2, with photometric redshift accuracy corresponding 
to the percent level, reaching 0.1% from spectroscopic measurements for a sub-sample of millions 
of them for 0.7 < z < 2.1 [34, 35]. The mapping of the galaxy clustering on such a volume, as well 
as the weak lensing shear associated to the induced ellipticity in galaxy images, represent Euclid 
primary targets in order to constrain the DE behaviour in the corresponding redshift interval. 

In order to be prepared to fully exploit the potential of the observations outlined above for the 
investigation of the dark cosmological components, an intense preparatory work aiming at a detailed 
simulation of the interplay between structure formation and CMB lensing is necessary. The starting 
point is the evaluation of the lensing pattern through the N-body simulations of structure formation. 
This has the potentiality to reproduce the full statistics of the lensing signal, relaxing the assumption 
of Gaussianity beyond the two point correlation function [36, 37, 38, 39]; in these works, the signal 
for the power spectrum, predicted semi- analytically, was correctly recovered, while all sky maps of 
deflection angle and projected lensing potential were constructed, as well as lensed CMB T, E and 
B anisotropics, down to an angular resolution of a few arcminutes. Similar studies, consisting in 
the simulation of the ISW through N-body LSS, have also been performed [40]. 

In this work we make a further progress along this direction. For the first time we conduct 
ray tracing studies through N-body simulations of non-ACDM cosmologies [for a recent review on 
N-body simulations in DE cosmologies, see 41]. The main focus of our analysis is the study and 
characterisation of the response of lensing to the underlying DE behaviour, focusing on the epoch 
at which the lensing is most effective. As representatives of various DE cosmologies, we adopt the 
ones implemented in the CoDECS^ publicly available suite of N-body simulations [42]. The latter 
include various realisations of dynamical DE models characterised by a direct interaction between 
DE and CDM particles. 

This work is organised as follows. In § 2 we describe the set of investigated cosmologies and 
their implementation through the CoDECS N-body simulations; in § 3 we review the basics of CMB 
lensing and the approach we use to trace CMB photons through the CoDECS structures. In § 4 we 
show our results. Finally, in § 5 we summarise our work and draw the concluding remarks. 



2 DE cosmologies in CoDECS 

As a first extension of CMB lensing studies to cosmological scenarios beyond ACDM, we focus 
on models where the role of DE is played by a classical scalar field (j) dynamically evolving in 
a self-interaction potential V {(/)), and characterised by a direct interaction with CDM particles. 
This class of cosmologies, generally referred to as "coupled DE" models (cDE hereafter) have been 
widely investigated in the literature as a possible way to alleviate the fine-tuning problems of the 
Cosmological Constant [see e.g. 43, 44, 45, 46, 47, 48, and references therein]. At the background 
level, cDE models are described by an interaction term between DE and CDM continuity equations, 
while baryons are kept uncoupled: 

Pe + 3i/pe= -y|/3c(0)g^, (2.2) 

Pb + 3Hpt = 0. (2.3) 

Here the coupling function /3c (0) sets the strength of the interaction and the sign of the quantity 
(/)/3c(0) defines the direction of the energy flow between the DE and CDM fluids. Although in the 
most general case the coupling function can directly depend on the scalar field itself, in the present 
work we will focus on a subset of the CoDECS cosmological models characterised by a constant 

^ www. euclid-ec.org 

^ www. mar cobaldi.it / CoDECS 
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coupling /3c{(f)) = /5c, and we will then omit the scalar field dependence of the coupling function in 
the remainder of the paper. In Eqs. (2.1)-(2.3) the overdot represents the derivative with respect 
to cosmic time t and Mp\ = l/y/SnG is the reduced Planck mass, with G the Newton's constant. 
According to Eqs. (2.1)-(2.2), different choices of the scalar field self-interaction potential V{<^) and 
coupling /3c determine a different evolution of the DE density = 0^/2 + and equation of 

state = 0^/2 — ¥{(/)) /p^^, which allows to put broad constraints on the model's parameters 

based on background observables only [see e.g. 43, 49, 50, 51, and references therein]. In particular, 
cDE models are in general characterised by a scaling regime in the matter dominated era during 
which DE provides a fraction of the total cosmic energy density, directly proportional to the square 
of the coupling, fi^ ^ 4/3^/3. This regime is known as the ^-Matter Dominated Epoch (or (/)-MDE, 
see Amendola [43]). The DE-CDM interaction also determines a time evolution of the CDM particle 
mass given by the integral of Eq. (2.2): 

^y^=-;|ft0, (2.4) 

where Mc is the mass of a CDM particle such that pc = Mcfic, with the number density ric being 
constant in comoving coordinates. Such a variation of the CDM particle mass induces in turn a shift 
of the matter-radiation equivalence redshift Zeq- All these effects allow to constrain the coupling 
strength, in particular, for the simplest class of models characterised by the constant couplings /3c 
considered in the present work. 

Besides providing a non-trivial evolution of the DE density and equation of state parameter, 
cDE models are also characterised by a long-range fifth-force, mediated by the DE scalar field and 
acting between CDM particles. In this respect, the relation between cDE cosmologies and modified 
gravity models, as e.g. scalar-tensor theories of gravity, has been extensively discussed by e.g. 
Pettorino & Baccigalupi [48]. The growth of linear density perturbations on sub-horizon scales can 
be described in the Fourier space and in the Newtonian gauge by the evolution equations for the 
CDM and baryon density contrast ^c,^: 



Sc = -2H 



Sc + ^7TG[pbSt + PcScTc] , (2.5) 



St = -2HSt + 4^G [ptSt + PcSc] , (2.6) 

where the factor Fc = (l + 4/3^/3) includes the effect of the fifth- force mediated by the DE scalar 
field for CDM density perturbations. The second term in the first square bracket at the right-hand 
side of Eq. (2.5), instead, represents the velocity dependent term arising for CDM perturbations as 
a consequence of momentum conservation, first described in detail in Baldi et al. [52]. Both these 
two additional terms contribute significantly to modify the evolution of density perturbations as 
compared to the ACDM cosmology and even to an uncoupled scalar field DE model. In particular, 
while the fifth-force is always attractive for the case - considered in the present work - of a single 
CDM particle species (while the same is in general not always true for multiple CDM families, see 
e.g. Brookfield et al. [53], Baldi [54, 55]), the friction term can act both as a proper friction or as a 
dragging effect, depending on the sign of the quantity (j)/3c. The interplay between these two terms 
can imprint very peculiar features on the growth rate of matter density perturbations, allowing to 
put tighter constraints on the coupling and potential functions and to possibly reconstruct their 
functional form. At the non-linear level, the acceleration equation of coupled (CDM) particles is 
characterised by the same two additional terms appearing in Eq. (2.5): 



n 



(2.7) 



with the only further complication that the particle's velocity and acceleration are no longer neces- 
sarily aligned, thereby significantly increasing the complexity of the possible interplay between the 
friction and fifth- force effects for the non-linear structure formation and evolution [see e.g. 56]. 
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The Codecs simulations [55], which represent the largest suite of cosmological N-body and 
hydrodynamical simulations of cDE models to date, have been obtained by making use of a modi- 
fication by Baldi et al. [52] of the widely used parallel N-body code GADGET [57] that allows to 
self-consistently simulate the evolution of structure formation processes in the context of generalised 
cDE cosmologies by including all the relevant features just discussed, from their specific background 
evolution to the mass variation of CDM particles, to the effects of the fifth-force and of the extra- 
friction acting on individual particles. Such modified code has been widely used in the past to 
investigate cDE cosmologies with both constant [52, 56, 58] and variable [59, 60, 61, 62] couphngs, 
as well as for different choices of the scalar field self-interaction potential. All the CoDECS simula- 
tions have been recently made publicly available through a dedicated web database^. At present, 
they include two distinct sets of runs, the L-CoDECS and the H-CoDECS. The L-CoDECS simulations 
consist of cosmological boxes of 1 comoving Gpc//i aside, filled with 1024^ CDM and 1024^ baryon 
particles. Both types of particles are treated with collisionless dynamics only, which means that 
baryonic particles are not considered as gas particles but just as a different family of collisionless 
particles distinguished from CDM. This is done in order to account for the effect of the uncoupled 
baryon fraction in cDE models which would not be correctly represented by CDM-only simulations. 
The mass resolution at z = for this set of simulations is Mc — 5.84 x IQ-*^^ Mq/Zi for CDM and 
Ml) — 1.17 X 10^^ M0//i for baryons, while the gravitational softening is set to = 20 comoving 
kpc//i, corresponding to 0.04 times the mean linear interparticle separation. The H-CoDECS simu- 
lations are instead adiabatic hydrodynamical simulations on much smaller scales, which we do not 
consider in the present work. 

The CoDECS suite presently includes six different cosmological scenarios: a fiducial ACDM 
cosmology taken as reference, three cDE models with a constant positive coupling /3c > and an 
exponential self-interaction potential: 

= Ae-^^ , (2.8) 

one cDE model with the same potential and an exponential coupling function: 

= /3oe^l^ (2.9) 

and one final cDE scenario with a constant negative coupling /3c < and a SUGRA [63] self- 
interaction potential: 

T/(0) = A(/)-^e^'/^ (2.10) 

The exponential potential of Eq. (2.8) is known to provide, in the absence of couplings, stable 
attractor solutions characterised by a constant ratio between the DE density and that of the dom- 
inant cosmological component [64]. When a coupling is active, however, the transient meta-stable 
scaling regime already discussed above and known as (/)-MDE occurs [43]. The SUGRA potential 
form derives instead from supersymmetric theories of gravity [65], determining a flattening of the 
potential at low redshifts in the case of no coupling; Baldi [55] has extensively investigated the 
same SUGRA model in the case of non-zero coupling [see also 66], showing that it determines a 
"bounce" of the DE equation of state it;^ on the cosmological constant barrier it;^ = — 1 at rela- 
tively recent times for initial conditions compatible with a static scalar field in the early Universe, 
^{z oc) = 0. Such cDE model, characterised by a constant negative coupling and by a SUGRA 
self-interaction potential, was then termed the "Bouncing cDE scenario", and has been shown to 
determine potentially observable effects on the number counts of massive clusters as a function of 
redshift [55]. All the models included in the CoDECS suite have the same background cosmological 
parameters at z = consistent with the latest results from WMAP-7 [67], summarised in Table 1, 
although the modified DE component makes the dynamics of both the cosmic background and the 
cosmological perturbations differ from ACDM. Linear density perturbations are normalised to the 
same amplitude at the redshift of the last scattering surface z\q ^ 1100, and evolved forward in time 
with the specific growth factor D^{z) of each model until the starting redshift of the simulations, 

^ www. marcobaldi.it / CoDECS 
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Parameter 



Value 



70.3 km Mpc 



-1 



0.226 
0.729 



2.42 X 10-*^ 



n. 



'S 



0.0451 
0.966 



Table 1. CoDECS cosmological parameters at z = 0: consistently with the notation adopted in the literature 
[67], Hq is the Hubble expansion rate at present, the density ratio with respect to the cosmological critical 
density of the species x; Ag and represent the amplitude and spectral index of the initial power law scalar 
perturbation spectrum, respectively. 

Zi = 99. The growth factor was obtained by numerically integrating Eqs. (2.5)-(2.6) along the re- 
spective background solution. Initial conditions of all the CoDECS simulations have been produced 
by setting up a random-phase realisation of the linear matter power spectrum obtained by running 
the public Boltzmann code CAMB^ for a ACDM cosmology with the parameters of Table 1, and 
rescaling particles displacements from a glass distribution using the different growth factors as de- 
scribed above. Therefore, all the simulations virtually share exactly the same initial conditions at 
the last scattering surface. 

As anticipated above, for the present work we restrict our investigation to a few selected 
models among the six cosmologies available within the CoDECS Project, with the aim of describing 
the main effects that cDE cosmologies generate on CMB lensing observables. Furthermore, the 
CoDECS suite will be soon extended to a wider range of cDE cosmologies, which will allow for 
an extensive application of our numerical approach in the near future. Our selected sample of 
models includes the reference ACDM scenario, the most strongly coupled model (/3c = 0.15) among 
those with a constant positive coupling plus an exponential potential, labelled "EXP003", and, 
finally, the Bouncing cDE model characterised by a constant negative coupling (/3c = —0.15) plus 
a SUGRA potential, labelled "SUGRA003". This choice is motivated by the fact that these two 
cDE models show very different background and linear perturbations evolutions, despite having the 
same absolute value of the coupling |/3c| = 0.15, and therefore the same strength of the associated 
fifth force, given by the factor Fc in Eq. (2.5) which for both models takes the value Fc = 1.03, 
corresponding to a 3% enhancement in the gravitational attraction for CDM particles. The different 
impact of these two models on observable quantities associated to the formation of linear and non- 
linear structures, including the CMB lensing signal, will then be entirely determined by the different 
dynamics of the DE scalar field and by the related evolution of the CDM particle mass variation 
and extra- friction term introduced in Eqs. (2.4) and (2.5), respectively. The viability of these 
models in terms of CMB observables has yet to be properly investigated, in particular for what 
concerns the impact of variable-coupling and bouncing cDE models on the large-scale power of 
CMB anisotropics. Although such analysis might possibly lead to tighter bounds on the coupling 
and on the potential functions than the ones allowed in the present work, here we are mainly 
interested in understanding the impact of cDE scenarios on the secondary anisotropics induced on 
CMB by large-scale structure lenses at late times, with a particular focus on the role played by 
non-linear effects, and we deliberately choose quite extreme values of the models parameters in 
order to maximize the effects under investigation. 

The relevant parameters of the cosmological models considered in this work are summarised 
in Table 2. For further details about the simulations we refer the interested reader to the CoDECS 
paper [42]. 

^ www . cosmologist .info 



6 



Model 



Potential 



a 



/3c 



Scalar 
field 
normalisation 



Wci){z 0) (j^{z 0) 



ACDM V{(j)) ^ A - - - -1.0 

EXP003 1/(0) = Ae-^^ 0.08 0.15 = 0) = -0.992 

SUGRA003 1/(0) = A0-^e^'/2 2.15 -0.15 0(z ^ oc) = -0.901 



0.809 
0.967 
0.806 



Table 2. Cosmological models from CoDECS which are considered in the present work and their specific 
parameters. 



3 CMB lensing 

CMB photons are deflected from an original direction on the last scattering surface to a direction n 
on the observed sky, so that the lensed CMB field is given by X(n) = -^(n^) in terms of the unlensed 
field X — T^Q^U, representing the total intensity and Stokes parameters for linear polarisation, 
respectively. As in previous works [36, 38], we consider here the case in which the change in 
the comoving separation of CMB photon trajectories induced by lensing is small compared to the 
comoving separation of the undeflected rays. In this case it is sufficient to calculate the relevant 
quantities for lensing, including the deflection angle, along the undeflected rays. This approach 
corresponds to the Born approximation. Adopting conformal time and comoving coordinates in a 
flat Friedmann Robertson Walker (FRW) geometry, the integral for the projected (along the line of 
sight) lensing-potential due to scalar perturbations with no anisotropic stress reads 



where r is the comoving distance, ~ 10^ Mpc is its value at the last-scattering surface, 770 is the 
present conformal time, $ is the gravitational potential generated by density perturbations, and 
[l/r]Vn is the two dimensional (2D) transverse derivative with respect to the line of sight pointing in 
the direction n = (??, (p). The vector is obtained from n by moving its end on the surface of a unit 
sphere by a distance |Vn^^(n)| along a geodesies in the direction of Vn^^(n). We assume |Vn^^(n)| 
to be constant between n and n', consistently with the Born approximation. The latter has been 
shown to hold down to the arcminute scale, taking into account the non-linearities corresponding 
to such scales [68, 69]. If the gravitational potential $ is Gaussian, the lensing potential is also 
Gaussian. However, the lensed CMB is non-Gaussian, as it is a second order cosmological effect 
produced by cosmological perturbations onto CMB anisotropics, yielding a finite correlation between 
different scales and thus non-Gaussianity. This is expected to be particularly important especially 
at small scales, i.e. tens of comoving Mpc or less, due to the the fact that non-linearities are already 
present in the underlying properties of lenses. It is typically more accurate to solve the integral in 
Eq. (3.2) directly to obtain the deflection angle instead of finite differencing the lensing potential 
[70, 71]. Semi- analytical techniques, such as LensPix and lenS2HAT [72, 73], exist for calculating 
Gaussian realisations of the lensing deflection angle, and therefore lensed CMB maps, starting from 
the lensing potential power spectrum, including non-linear corrections on small scales [74]. This 
approach has been used to study and validate CMB lensing templates obtained via ray tracing 
through the Millennium Run N-body simulation at the level of the two point correlation function 
[38], finding significant deviations from the semi-analytical expectations at multipoles of order 10^. 
On large scales, corresponding to about 1/3 of a degree in the sky or more, due to the limited box- 
side size, L = 500 Mpc/h comoving, of the Millennium Run, the lensing signal was not correctly 




(3.1) 



while the corresponding deflection angle integral is 




(3.2) 
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reproduced. Therefore, Carbone et al. [36, 38], as a feasibility study, imported synthetic maps 
of lensing potential, deflection angle and lensed CMB generated through Lenspix, filtered on large 
scales, in order to achieve the complete angular power spectrum (for further details see Carbone et 
al. [36, 38]). Even if the CoDECS simulations still represent only a portion of the Hubble volume, 
this technique is not applied in the present analysis since we are interested to the pure N-body 
signal and its response to the underlying DE behaviour. In the next two sub-sections we outline 
our ray-tracing and sky projection procedures. 

3.1 Ray tracing 

The lensing potential maps obtained through ray-tracing across the CoDECS structures have been 
generated as described below. 

First, we pre-compute and store the gravitational potential grids from the simulation snapshots, 
taking care of achieving a dense sampling in redshifts. In Fig. 1, we show the grid spacing x 
representing the comoving distances of the snapshots as a function of the redshift, together with 
their ratio to the CoDECS box-side size, Lbox = 1 Gpc/h. Notice that we had to re-sample the 
simulation outputs with respect to the public version, in order to achieve a proper sampling in z for 
CMB lensing analyses. The original redshift sampling of the CoDECS snapshots is shown in red in 
the upper and lower panels of Fig. 1. 

Second, in order to produce mock lensing potential maps that cover the past light cone over the 
full sky, we stack the gravitational pre-computed potential grids around the observer following the 
technique developed in [36]. 

Third, the total volume around the observer up to is divided into spherical shells, each one with a 
thickness corresponding to Lbox- All the boxes falling into the same shell are translated and rotated 
with the same random vectors generating a homogeneous coordinate transformation throughout the 
shell, while the randomisation vectors change from shell to shell [36]. 

Finally, the peculiar gravitational potential at each point along a ray with direction n is spatially 
interpolated from the CoDECS grid which possesses a spatial resolution of Lbox/2560 ^ 400 kpc//i. 
The deflection angle is computed along the line of sight as well, by numerically evaluating the 
gravitational potential gradient and interpolating at each point along the line of sight. By varying 
the direction of the light-ray along which the integration is performed, we construct all sky maps 
of the lensing potential. 

3.2 Sky projection 

We adopt the Hierarchical Equalised Altitute Pixel (HEALpix^^) scheme for the pixelisation of 
the sphere [75], and choose a resolution of about 1.72 arcminutes, corresponding to a resolution 
parameter A^side — 2048; in an HEALPix environment, the latter choice corresponds to angular 
multipoles of at least Imax — 2 • A^nside = 4106, before pixel effects start to kick in. The chosen 
angular resolution has been checked as follows. In the left panel of Fig. 2 we compare the map 
angular resolution of 1.72' (solid black line) with the effective angular resolution 6(z) corresponding 
to the intrinsic grid spacing 400 kpc//i) of the three dimensional gravitational potential field as 
a function of redshift. We observe that at redshifts z < 0.3 the chosen angular resolution of 1.72' 
becomes smaller than 0{z). However, since most of the lensing signal comes from higher redshifts, 
the latter effect is negligible [36]. On the right panel of Fig. 2, we make a similar comparison, 
considering the maximum multipole, lmax{z) — 180/0(2:), that can be exploited given the intrinsic 
grid resolution. For most of the redshifts, the latter is larger than Imax — 4106. Nonetheless, 
we adopt the conservative choice of limiting our analysis to / < 3000 where Poisson noise from 
low-redshift potential integration starts to dominate the simulated CMB lensing signals. 

4 Results 

We now proceed in the illustration of our results. In Figs. 3, 4, and 5 we show the two components 
of the deflection angle field a along the 9 and (j) HEALPix spherical coordinates (azimuth and 

^"^http:/ /healpix.jpl.nasa.gov 
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Figure 1. Top: comoving distance of the gravitational potential snapshots extracted from the CoDECS 
simulations. Bottom: ratio of the comoving spacing with respect to the box side. The red data correspond 
to the original spacing in CoDECS, which was re-computed to become the one corresponding to the green 
points, in order to implement the CMB ray-tracing across the gravitational potential of the simulations. 
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Figure 2. Left: comparison between the pixel size = 1.72' of the lensing-potential sky maps (dotted 
black line) and the redshift dependent angular resolution (solid blue line) corresponding to the intrinsic 
spacing (~ 400Mpc//i) of the three-dimensional gravitational potential grids extracted from CoDECS. Right: 
comparison between the maximum multipole Imax corresponding to 2 • A^side and the redshift dependent 
maximum multipole Imaxi^) = 180/^ corresponding to the resolution of the potential grids. 



elevation angles, respectively), its magnitude \a\, as well as the projected lensing potential, for the 
three models described in §2, namely ACDM, EXP003, SUGRA003, respectively. For these 
models, the pattern across the sky is similar in morphology, reflecting the same initial conditions in 
the CoDECS simulations. The amplitude of the lensing signal, though, depends on the integral over 
the subsequent structure formation epoch, weighted with the lensing cross section, (r* — r)/(r*r), 
in Eq. (3.1). Already at the level of the maps, the increase in the CMB weak lensing power for 
the EXP003 model can be seen when looking at the maximum value assumed by \a\ in the lower 
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LCDM: total deflection— angle map along the theta— direction 



LCDM: total deflection— angle map along the phi— direction 
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Figure 3. Maps of the relevant quantities for CMB weak-lensing for the ACDM model: deflection-angle 
components along the azimuth and elevation angles (top panel: left and right, respectively), its magnitude 
(in radians) and the dimensionless lensing potential (bottom panel: left and right, respectively). 



EXP003: total deflection— angle map along the theta— direction EXP003: total deflection— angle map along the phi— direction 
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Figure 4. As in Fig. 3 for the EXP003 model. 
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SUGRA003: total deflection— angle map along the th eta— direction 



SUGR4003: total deflection— angle map along the phi— direction 




SUGRA003: total deflection— angle modulus map SUGR4003: total lensing— potential map 




Figure 5. As in Figs. 3-4 for the SUGRA003 model. 



left panel of Fig. 4, as compared to the same quantity in Fig. 3, and reflects the behaviour of 
structure formation, affected by the different growth factor and density perturbation amplitude 
with respect to the ACDM case. On the other hand, the amplitude of the lensing effects for the case 
of SUGRA003 becomes closer again to the ACDM model, although as we discuss below, a finite 
and macroscopic difference remains also in this case, reflecting the markedly different background 
and dynamic evolution at the relevant epochs for CMB lensing. 

More quantitatively, the mean value of the deflection angle modulus in the simulated maps 
is given by mean(|a|) = 1.75^ 2.12', 1.66^ and the standard deviation is equal to cr(|Q:|) = 
0.92', 1.11', 0.87', for the adopted models ACDM, EXP003 and SUGRA003, respectively. For com- 
parison, a synthetic Gaussian map, produced with the lensing potential power spectrum generated 
with CAMB (i.e. taking into account also the scales not covered by the simulation box size) using 
a WMAP-7 ACDM cosmology [67], is characterised by cr(|a|) = 1.2'; therefore the lack of scales 
larger than the simulation box results in a suppression in cr( |a|) of about 23% in the simulated CMB 
lensing signal. 

As explained in §2, in the CoDECS simulations, baryon particles are assumed to be collisionless. 
Nonetheless, we can safely neglect the contribution of gas dynamics and feedback on the scales 
considered in this work, since they are limited by the intrinsic resolution ^ 400 kpc//i of the 
gravitational potential grids, while, as shown in [76], baryonic physics alters the matter power 
spectrum on much smaller scales. 

Therefore, in Figs. 6-8 we show collisionless baryon contributions to the quantities relevant 
for CMB lensing. We observe that baryons mostly follow the dynamics of CDM particles, which 
are the dominant non-relativistic component. For this reason, the morphology of the lensing po- 
tential and deflection angle maps is markedly similar between baryons and CDM, and the general 
trend, described above for the total matter contribution, is reproduced individually by each of the 
components. On the other hand, the contributions to the signal amplitude depend on the relative 
abundances. In particular, for ACDM, EXP003 and SUGRA003, respectively, baryons contribute 
to a by mean(|a|) = 0.29', 0.34', 0.30', and a{\a\) = 0.15', 0.18', 0.16', while CDM particles give 
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LCDM: deflection— angle map along the theta— direction for coUisionless baryons 



LCDM: deflection— angle map along the theta— direction for CDM 




LCDM: deflection— angle map along the phi— direction for coUisionless baryons 
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LCDM: deflection— angle map along the phi— direction for CDM 
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LCDM: deflection— angle modulus map for collisionless baryons 



LCDM: len sing— potential map of collisionless baryons 
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LCDM: deflection— angle modulus map for CDM 




LCDM: lensing— potential map of CDM 




Figure 6. Maps of the relevant quantities for CMB weak-lensing in the ACDM case: from top to bottom, 
deflection angle components along the azimuth and elevation angles, magnitude and lensing potential (in 
radians and dimensionless, respectively), for baryons (left) and CDM (right). 



mean(|Q:|) = 1.46^1.78^1.36^ and cr(|a|) = 0.76^ 0.93^ 0.71^ Therefore, baryons produce only 
^ 1/3 of the total rms of the deflection angle modulus. This can be seen, for example, by looking 
at Fig. 9: here we plot the lensing potential produced by both baryons and CDM components, 
for the three models discussed above in the text. First, we see that there is a difference in power 
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EXP003: deflection— angle map along the theia— direction for coUisionless baryons 



EXP003: deflection— angle map along the theta— direction for CDM 




of order ^ 0{1) between the CDM and the baryon spectra. Moreover, the trends of the EXP003 
lensing potential C/, which are higher than the corresponding spectra in the ACDM model at all 
scales, both for the CDM and baryon cases (dashed and dot-dashed hnes, respectively, in Fig. 9), 
are well expected, given the higher matter power spectrum normalisation at the present epoch with 
respect to the ACDM cosmology. 
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SUGRA003: deflection— angle map along the tetha— direction for coUisionless baryons 



SUGRA003: deflection— angle map along the tetha— direction for CDM 




On the contrary, for the SUGRA003 model, baryons and CDM have a different and unexpected 
trend with respect to the ACDM and EXP003 ones. Indeed, for SUGRA cosmologies, the CDM 
lensing potential power (3-dot-dashed hne) hes below the ACDM one (solid hne), while the baryon 
lensing potential power (long-dashed line) is slightly above the ACDM one (dotted line). 
In all the three cases, the differences among the baryon lensing potential power spectra are visibly 



-14- 



10 



-6 




10 



-10 



baryons-SUGRA003 



baryons-LCDM 
baryons-EXP003 



10 



-11 



10 



100 

1 



1000 



Figure 9. Lensing potential power spectra of the baryon and CDM components, for the models considered 
in the text. 



less marked than for the CDM components. Therefore, even if collisionless, baryon dynamics does 
not perfectly follow the CDM dynamics. The different behaviours between the baryon and CDM 
cases, in the three cosmologies considered, are due to the DE coupling with CDM, which is absent 
for baryons. In addition, at the end of this Section, we give a physical explanation to the unexpected 
CDM lensing potential trend for the SUGRA003 case. 

In the upper panel of Fig. 10 we show the total lensing potential power spectra extracted from 
the simulations, for the ACDM, EXP003 and SUGRA003 models. Several features can be noticed 
in this figure; we start with comments related to the properties of the simulations, regardless of 
the cosmological model, and then we describe the different features reflecting the underlying cDE 
models. 

First, the scale in which CoDECS allow to trace real structures in the lensing pattern is broader 
by a factor of ^ 2 with respect to the results of [36, 38], obtained via ray-tracing across the 
Millennium Run simulation. In particular, we have checked that, in the ACDM case, CoDECS 
reproduce correctly the lensing potential power on scales / > 200. This is due to its box side, which 
is doubled with respect to the Millennium Run, able to reproduce the lensing signal only for / > 400. 
Second, the plot shows how the lensing power faithfully traces differences in structure formation 
at the epoch in which the lensing cross section for CMB lensing is maximum, i.e. z ^ 1 zb 0.5, as 
anticipated in earlier works [27]. This effect is described by the different amplitudes and shapes of 
the lensing potential power spectra as represented by the solid black line, the dashed blue line and 
the dot-dashed red hne, for the ACDM, EXP003 and SUGRA003 models, respectively. 

The upper panel of Fig. 10 has to be compared with Fig. 11, in which we show the 
semi-analytical predictions obtained with the implementation of the publicly available Code for 
Anisotropics in the Microwave Background (CAMB^^), adapted to interacting dark energy cos- 
mologies [51]. 

From this comparison, we notice that, on the scales probed by the simulations, / > 200, there is 
a qualitative good agreement between semi-analytical results and simulated data, but, except for 
the ACDM case, already at multipoles / > 400, where the non-linear structure evolution starts to 

^^See http://camb.info 



15 




2 
Q 

CO 

o 
u 

Q 
O 
O 



L5 - 



LO 



0.5 - 



0.0 



CDM + baryons EXP003 
CDM+baryons SUGRA003 
CDM+baryons EXP003 normalised to LCDM 



10 



100 



1000 



Figure 10. Upper panel: lensing potential angular power spectra for the models considered in the present 
work. Lower panel: ratios of the simulated EXP003 and SUGRA003 lensing potential power spectra with 
respect to the ACDM one. The EXP003 signal is also shown after re-scaling to the matter power normalisation 
of the ACDM model at present, as explained in the text. 



affect the lensing signal, the semi-analytical curves deviate from the simulated ones, due to the 
Halofit [74] non-linear implementation present in CAMB, which has been tested against N-body 
simulations only in ACDM cosmologies. 

We avoid to compare in the same plot the semi-analytical expectations with the simulated results 
from Codecs. In fact, the integration of background equations for coupled quintessence models 
requires an iterative routine that allows to tune the initial conditions such that the wished back- 
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Figure 11. Semi- analytic expectations of the lensing potential angular power spectra for the ACDM, 
EXP003, and SUGRA003 models, obtained with a modified CAMB version, as described in the text. 



ground cosmological parameters at present are recovered. For the SUGRA models, the presence of 
the bounce increases the level of fine tuning and makes the background evolution more unstable 
with respect to small variations of the coupling and potential parameter. As a consequence, the 
background values obtained in the SUGRA003 model for CAMB are not exactly the same as the 
ones used for the CoDECS background. For this reason, in order to compare the semi-analytical 
expectations of the lensing potential among the three cosmological models considered in this work, 
and check if they confirm qualitatively the Ci trends extracted directly from CoDECS, we have de- 
cided to produce the corresponding ACDM and EXP003 spectra using exactly the same background 
values chosen in CAMB for SUGRA003, instead of adopting the CODECS background values. 

In the lower panel of Fig. 10, we show the ratios of the lensing potential spectra with respect to 
the ACDM signal, as extracted from the simulations, and, for the EXP003 case, we also display the 
same ratio normalised to cr8(ACDM)5 obtained dividing the EXP003 lensing potential Ci by ^l^^^xprn'^) 
and multiplying it by (^\i^^Qj>,y\^^ (see Table 2). In particular, this model shows a power excess on all 
scales / > 200 of ^ 50% with respect to the ACDM spectrum (dashed blue line in the lower panel 
of Fig. 10). As we can observe from the dotted green line in the same plot, this excess of power is 
mostly due to the different matter power spectrum normalisation at z = between the ACDM and 
EXP003 cosmologies arising as a consequence of the enhanced linear growth. Nonetheless, on non- 
linear scales, / > 1000, a residual excess is present, produced by the different non-linear dynamics 
and evolution of the CDM perturbations, given the matter coupling with the DE. At / ^ 3000 
non-linearities contribute to the power excess in the EXP003 lensing signal by ^ 10%. This effect 
can not be reproduced with the use of semi-analytical techniques in cDE models, being the product 
of the LSS non-linear evolution accurately described only via N-body simulations. Furthermore, it 
would be misleading to conclude that the as degeneracy between ACDM and EXP003 breaks only 
beyond the arcminute scale, i.e. where the non-linear structure evolution is more efficient, especially 
due to the coupling between DE and CDM [54]. This is what one could wrongly infer owing to 
the limited box size of the adopted N-body simulations. In fact, having a look at the semi-analytic 
predictions in Fig. 11, which reproduce the lensing signal on scales / < 200, it is easy to observe that 
the as degeneracy between the EXP003 and ACDM cosmologies is completely broken if also scales 
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Figure 12. The redshift evolution of the three terms entering the expression for the lensing power spectrum 
amphtude, and their combination, as compared to the ACDM case. As one can see from the plots, the 
different terms have very different evolutions that partly compensate with each other, resulting in an overall 
evolution of the lensing power amplitude as compared to ACDM displayed in the lower-right panel. The 
latter clearly shows how an enhanced lensing signal is expected for the EXP003 model, while a moderately 
suppressed lensing power is predicted for SUGRA003, despite for this model the amplitude of matter density 
perturbations is always larger than in ACDM. 



larger than the simulation box are taken into account. This feature will be transferred, after lensing, 
to the first CMB temperature and polarisation peaks, allowing to clearly distinguish between the 
two models even after normalisation to the same a^. Lensed T, E, and B CMB spectra in cDE 
cosmologies will be accurately treated in a future work. 

Finally, as the dot-dashed red lines in the upper and lower panels of Fig. 10 clearly show, 
and as anticipated above, the lensing effect produced by the SUGRA003 model shows a completely 
opposite trend with respect to the EXP003 one. In fact, in this case, the combination of modified 
geometry and abundance of lenses determines a loss of power, remarkably at the 10% level on all 
the scales spanned by the simulations. We have also checked that, as a consequence of SUGRA003 
being similar to ACDM in perturbation power at present, this effect stays unchanged after as 
normalisation. Such phenomenology might appear counter-intuitive as the SUGRA003 model is 
expected to always have a higher amplitude of matter density perturbations as compared to ACDM, 
except for z = ^cmb and z = [42]. One would intuitively expect also a larger lensing potential 
arising as the integral over the line-of-sight of a density power spectrum with a higher amplitude 
than the standard ACDM cosmology. However, this is not the case due to the interplay between 
the evolution of background quantities and density perturbations in determining the lensing power 
spectrum P^. In fact, within the small-scale Limber approximation, the latter is related to the 
matter density power spectrum through the relation: 

= (4.1) 

such that the excess with respect to the standard ACDM cosmology can be written as the combi- 



- 18 - 



nation of three main factors: 
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(4.2) 



where is the hnear growth factor normahsed at recombination, and where we have taken into 

account that ah the cosmological models under investigation share the same hnear matter power 
spectrum at last scattering. 

These three different terms are shown in Fig. 12, together with their combined effect on the 
lensing power spectrum amplitude. While the evolution of the growth factor alone would suggest 
a higher lensing efficiency due to the larger amplitude of density perturbations, the remaining 
two terms compensate for this effect in a non-trivial way, determining an overall evolution of the 
expected lensing power amplitude fully consistent with the results obtained through our ray-tracing 
procedure. Clearly, such estimate is expected to fail at highly nonlinear scales, where in fact our 
fully nonlinear treatment is capable to capture additional features that are not accurately predicted 
by the ratio of Eq. (4.2). 

5 Conclusions 

The modification of the anisotropics in the CMB induced by forming cosmological structures through 
gravitational lensing is one of the most important topics in modern cosmology. The lensing signal 
probes the expansion rate, and thus the abundances of the dark cosmological components, includ- 
ing DE, through the distribution of lenses at an epoch which corresponds to the onset of cosmic 
acceleration. 

In view of the existing plans for ground and space based wide field and deep observations of LSS, 
large N-body simulations are being produced or planned, involving DE models as well as scenarios in 
which gravity is modified, in order to simulate with high accuracy the cosmic acceleration in theories 
alternative to the cosmological constant scenario. CMB anisotropy data, on the other hand, are 
available over the entire sky, with angular resolution reaching a few arcminutes, and sensitivity of 
a few micro-Kelvin, in total intensity and polarisation. 

In this context, it is important to produce simulations of the expected CMB lensing pattern 
in the aforementioned scenarios, for two reasons. A first one concerns testing and developing 
our capability of implementing N-body based simulations of CMB lensing, involving ray tracing 
through simulated structures. A second one consists in studying and quantifying the sensitivity of 
CMB lensing variables to the underlying cosmological model, concerning the modelling of DE, in 
particular. 

In this work, we progressed on both aspects, by simulating the CMB weak lensing potential 
through N-body simulations featuring various DE models characterised by modified expansion his- 
tories, modified cosmological perturbation dynamics, and the presence of a coupling between CDM 
particles and DE, in comparison to the expectations of the ACDM model. We constructed suitable 
grids of the gravitational potential output from N-body simulations, and implement ray tracing in 
the Born approximation. We derived maps of the deflection angle and lensing potential, verifying 
our capability of tracing N-body structures in the CMB lensing signal on a range of scales going 
from the arcminute to the degree scale, dictated, as expected, entirely by the simulation resolution 
and extension. Following well known geometrical properties of CMB lensing, which has its cross 
section maximum at redshifts of about z 2^ 1 ± 0.5, we studied how the simulated lensing signal 
traces differences in the background and perturbation dynamics at the corresponding epoch. We 
considered two DE models, EXP003 and SUGRA003, the first featuring an increased structure for- 
mation rate with respect to the ACDM at all epochs, and the second characterised by a matter 
perturbation normalisation close to the ACDM at the initial and present epochs, being different 
only at intermediate redshifts. We show how the lensing variables faithfully record the differences 
of the models with respect to the ACDM scenario at the onset of acceleration, in particular for the 
second scenario, degenerate with the ACDM at present and at the initial time. Correspondingly 
to the different amplitudes and dynamics of the distribution of lenses in the various models, for 
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the SUGRA003 cosmology the lensing potential power results to be of order 10% smaller than the 
corresponding ACDM signal on all the scales covered by the simulation. This is a quite unexpected 
result, and it is due to the compensation between the larger power of the matter perturbations 
and the suppressing effect of the background and matter density evolution characterising this spe- 
cific class of cDE models. The CoDECS simulations store the gravitational potential of baryons 
and CDM particles separately; therefore, we inspected the results for the baryon component only, 
finding a substantial compensation as for the CDM particles, with the power being slightly higher 
than ACDM in this case. We attributed this occurrence to the coupling with DE, which concerns 
CDM particles, and not baryons. For the EXP003 model, we find an opposite result, i.e. we find an 
excess of order 10%, with respect to the ACDM case, for the lensing potential power on multipoles 
/ > 1000, after normalisation of the signals to the same as- On the large scales, i.e. on multipoles 
/ < 20 instead, the semi-analytical expectations predict that the signals tend to converge to the 
same amplitude independently of the underlying cosmology, owing to the same curvature pertur- 
bation power amplitude assumed for the initial conditions of the N-body simulations. This means 
that the degeneracy between the ACDM and EXP003 model can be broken if a large interval of 
scales is analysed in the observations. 

These results confirm the relevance of CMB lensing as a probe for DE at the early stages 
of cosmic acceleration, and demonstrate the reliability of CMB lensing based on the existing N- 
body simulations, in the context of DE studies. At the same time, our study clearly indicates 
paths forward, especially in terms of box size and alternative simulated cosmologies, in order to 
be able to access super-degree angular scales, as well as accurately mapping the lensing pattern on 
arcminute and sub-arcminute angular scales; the latter in particular represents a crucial regime for 
CMB lensing, as it gets to dominate with respect to the decaying CMB anisotropics of primordial 
origin due to diffuse damping at last scattering, and at the same time involving non-linear physics 
where semi- analytic modelling may be either unavailable or inaccurate. In this respect, N-body 
simulations of actual structure formation may be the only way to well characterise DE effects in 
such a regime. 
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